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with q colors, and define W = lim^—^ P(G, q) 1 ^. Besides their mathematical 
interest, these functions are important in statistical physics. We give a com- 
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families of graphs, including strips of regular lattices with various boundary 
conditions and homeomorphic expansions thereof. Generalizing to q £ C, we 
determine the accumulation sets of the chromatic zeros constituting the con- 
tinuous loci of points on which W is nonanalytic. Various families of graphs 
with the property that the chromatic zeros and/or their accumulation sets (i) 
include support for Re(q) < 0; (ii) bound regions and pass through q = 0; 
and (iii) are noncompact are discussed, and the role of boundary conditions is 
analyzed. Some corresponding results are presented for Potts model partition 
functions for nonzero temperature, equivalent to the full Tutte polynomials 
for various families of graphs. 
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I. INTRODUCTION 



We consider connected graphs G without loops or multiple bonds and denote the number 
of vertices as n = v(G), the edges as e(G), the girth as g, the number of circuits with this 
minimal length as k g , and the chromatic number as The chromatic polynomial P(G, q) 

counts the number of ways that one can color the graph G with q colors such that no two 
adjacent vertices have the same color [l[ (for reviews, see [[|- The minimum number 
of colors needed for this coloring, i.e., the smallest positive integer q for which P(G,q) is 
nonzero, is the chromatic number, x{G)- Besides its intrinsic mathematical interest, the 
chromatic polynomial has an important connection with statistical mechanics since it is the 
zero-temperature partition function of the g-state Potts antiferromagnet (AF) ||[7j on G: 

P(G,q) = Z(G,q,T = 0) PAF (1.1) 

We shall consider recursively defined families of graphs here, i.e., roughly speaking, those 
that can be formed by successive additions of a given subgraph. A precise definition of a 
recursive family G m depending on a positive integer parameter m is a sequence of graphs for 
which the chromatic polynomials P(G m , q) are related by a linear homogeneous recurrence 
relation in which the coefficients are polynomials in q |],^|. These include strips of regular 
lattices with various boundary conditions, chains of polygons linked in various ways, joins 
of such graphs with another graph such as a complete graph K p , and families obtained from 
these by modifications such as removal of some edges or additions of degree-two vertices, 
i.e., homeomorphic expansions. These families of graphs may depend on several parameters 
(e.g., width, length, number of homeomorphic expansions, etc.), and there can be several 
ways in which one can obtain the limit n — > oo. Let us concentrate on one such parameter, 
such as the length of a strip of a regular lattice of fixed width or polygon chain, m, so that 
n is a linear function of m. 

Just as the chromatic polynomial counts the total number of ways of coloring the graph 
G with q colors subject to the above constraint, so also it is useful to define a quantity that 
measures the number of ways of performing this coloring per site, in the limit where the 
number of vertices goes to infinity. Denoting the formal limit {G} = lim^oo G, we havef] 



At certain special points q s (typically q s = 0,1, ..,x{G)), one has the noncommutativity of limits 
lim g ^ gs lim n _,oa P(G, q) 1 ^ 71 ^ \im n - >00 lim q ^ qs P(G, q) 1 ^ 1 , and hence it is necessary to specify the order 
of the limits in the definition of VF({G}, q s ) ||. We use the first order of limits here; this has the advantage 
of removing certain isolated limits here; this has the advantage of removing certain isolated discontinuities 
in W. 
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W({G},q)= limP(G,g) 1 / n (1.2) 

n—*oo 

With the order of limits specified in the footnote, this limit exists for the recursive families 
of graphs considered here, as will be discussed further below. The function W is the ground 
state degeneracy per vertex in the n — > oo limit. The quantity So = ks In W, where ks is 
the Boltzmann constant, is defined as the ground state entropy. The Potts antiferro magnet 
has the interesting feature that it exhibits nonzero ground-state entropy So ^ (without 
frustration) for sufficiently large q on a given lattice or graph and is thus an exception to 
the third law of thermodynamics |1(J. This is equivalent to a ground state degeneracy per 
site W > 1. For example, for the square lattice, W = (4/3) 3/2 0. We recall that, with 



n = v(G), a general form for the chromatic polynomial of a connected graph G is 



ra-1 



p(g, q ) = Tr,(-iYK-j<r' (i.3) 

where h n -j are positive integers, with [[4|,|T2"f h n _j — ( e J for < j < g — 1 (whence h n = 1 
and /i n _i = e) and fr„-( g -i) = ( g f J - fc s - 

Although in the Hamiltonian formulation of the g-state Potts model or the mathematical 
context of coloring a graph with q colors, q must be a non-negative integer, once one has the 
function P(G, q), one can generalize the quantity q from Z + to C. A subset of the zeros of 
P in the q plane (chromatic zeros) may form a continuous accumulation set in the n — » oo 
limit, denoted B, which is the continuous locus of points where W({G},q) is nonanalytic 
(B may be null, and W may also be nonanalytic at certain discrete points). The maximal 
region in the complex q plane to which one can analytically continue the function ^({G}, q) 
from the range of positive integer values where there is nonzero ground state entropy is 
denoted R\. The maximum value of q where B intersects the (positive) real axis is labelled 
q c ({G}). Thus, region Ri includes the positive real axis for q > q c ({G}). We have calculated 
B for a variety of families of graphs and shall present rigorous results, some observations, 
and some conjectures here. Part of the interest in this area of research is that it combines, 
in a fruitful way, both theoretical physics and three areas of mathematics: graph theory, 
complex analysis, and algebraic geometry. In addition to the works already cited, some 
previous relevant papers are [13|- pl| . 



The recursive families of graphs considered here include strips of regular lattices with 
various boundary conditions, chains of polygons linked in various ways, joins of such graphs 
with another graph such as a complete graph K p , and families obtained from these by modifi- 
cations such as removal of some edges or additions of degree- two vertices, i.e., homeomorphic 
expansions. These families of graphs may depend on several parameters, and there can be 
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several ways in which one can obtain the limit n — > oo. Concentrating on one such parame- 
ter, such as the length of a strip of a regular lattice or polygon chain, m, and denoting the 
graph as (G s ) m , a general form for P((G s ) m ,q) is 

N a 

P((G s ) m , q ) = J2cM*M m (1-4) 

i=i 

where Cj(q) and the N a terms aj(q) (which we shall also denote equivalently as \j(q)) depend 
on the type of strip graph G s but are independent of m. We define a term ae as "leading" (£) 
if it dominates the n — > oo limit of P(G, q). The locus E occurs where, as one changes q, there 
is an abrupt, nonanalytic change in W as the leading terms ae in eq. (|L4|) changes. Hence, E 
is the solution to the equation of degeneracy of magnitudes of leading terms, |o^| = |<v|. It 
follows that |W| is finite and continuous, although nonanalytic, across E. In general, for the 
families of graphs considered here, E is an algebraic curve, since it is the locus of solutions to 
the equation \ae\ = \ay\ and the terms ae, ay are algebraic functions of q. From (|1.4j ), with 
the ordering as given in the footnote, one can show constructively that the limit (|1.2| ) exists; 
it is given by W({G s },q) = (ae) 1 for q G and q)\ = |c^|* for q £ Ri, where 

linim^oo m/n = t and ae is the term among the cij, j = 1, ...,N a with maximum magnitude 
at the given value of q. Thus, a number of interesting questions can be asked and answered 
about its properties. 

Concerning the coefficients Cj, we note that 

Ci(q) = 1 if ae{q) is dominant in Ri (1-5) 

This is proved by using the fact that region R\ is defined to include the region of large 
real q, and then requiring the expression ( |1.4|) to match the general expression ( |1.3| ) and 
satisfy the property that the coefficient of the leading term, q n , is unity. For some families of 
graphs, the Cj and aj are polynomials in q. However, there are also families of graphs, such 



as the L y = 3 cyclic and Mobius strips of the square and kagome lattices |36| , |39|| (see below 
for notation) which have nonpolynomial algebraic a^'s with polynomial c/s, and families of 
graphs, such as certain strip graphs of regular lattices with free transverse and longitudinal 



boundary conditions pUlBTJ, as well as the L y = 2 Mobius strips of the triangular lattice and 



asymmetric homeomorphic expansions of the L y = 2 Mobius square strip |36| , p9| for which 
some Oj's and their respective Cj's are both algebraic, nonpolynomial functions of q. When 
the Oj's are algebraic nonpolynomial functions of q, this happens because these are roots 
of an algebraic equation of degree 2 or higher. There are then two possibilities: first, the 
roots may enter in a symmetric manner, and, because of a theorem on symmetric polynomial 
functions discussed below, their coefficients are all equal and are polynomial functions of q. 
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In particular, if one of these a/s is leading in Ri, then the coefficient functions for all of the 
roots are equal to unity, in accordance with ( |1.5|) . This happens, for example, for the L y = 3 
cyclic and Mobius strips of the square lattice, where the leading term, a sq ^ = \ sq fi in eq. 
( |4. 1|) is the root of a quadratic equation (see eq. Q4.2Q ) and enters the chromatic polynomial 
together with the other root of this equation, each with the same coefficient, which is unity. 
A second possibility is that the nonpolynomial algebraic roots enter in (|1.4p in a manner 
that is not a symmetric function of these roots; in these cases, the corresponding coefficient 
functions are also nonpolynomial algebraic functions of q. For example, in many strips with 



free transverse and longitudinal boundary conditions f3(J.|3T| and in the the Mobius strip of 
the triangular lattice and the homeomorphic expansion of the Mobius strip of the square 
lattice, one finds |36j that two of the Oj's are roots of a quadratic equation of the form 



a± = R ± y/~S, and they occur in the chromatic polynomial in the form (a+) m — (a_) m ; for 
these families, this combination is multiplied by the factor l/(a+ — a_), i.e., the coefficient 
functions contain the respective factors ±l/Vo. In all cases, of course, these combine to 
yield a P(G, q) that is a polynomial in q. 

Two other general result that follow from eqs. (|1.2| ), ( |1.3|) , and ( |1.4|) are as follows. As 
m, n — » oo, let r = n/m, which depends on G s . Then 

ae ~ q r for q £ R\ , \q\ ^ oo (1.6) 

W = (a e )r for q e Ri (1.7) 
so that limigi^oo; geRl q~ x W = 1. 

II. CALCULATION OF CHROMATIC POLYNOMIALS 

For a recursively defined family of graph G s our calculational method is to use the 
deletion-contraction theorem iteratively. For some families of graphs we have directly solved 
for the chromatic polynomials from the recursion relations that follow from this iterative pro- 
cedure. A related method that that we have found to be useful employs a generating function 



and is described below p0| , |31| . Let us consider strips of various lattices with arbitrary length 
L x = m vertices and fixed width L y vertices (with the longitudinal and transverse directions 
taken to be x and y). The chromatic polynomials for the cyclic and Mobius strip graphs of 
the square lattice were calculated for L y = 2 in by means of iterative deletion-contraction 
operations, and subsequently via a transfer matrix method in JT^] and a coloring matrix 



method in [0 (see also ||40fl). We have extended this to the corresponding cyclic |3(J and 
Mobius |3£| L y = 3 strips. For each graph family, once one has calculated the chromatic 
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polynomial for arbitrary length, one can take the limit of infinite length and determine the 
accumulation set B. After studies of the chromatic zeros for L„ 



2 in |,22,23 , W and 



B were determined for this case in H and for L,, 



3 in 



An interesting question 



concerns the effect of boundary conditions (BC's), and hence graph topology, on P, W, and 
B. We use the symbols FBC y and PBC y for free and periodic transverse boundary conditions 
and FBCa;, PBC^, and TPBCa, for free, periodic, and twisted periodic longitudinal boundary 
conditions. The term "twisted" means that the longitudinal ends of the strip are identified 
with reversed orientation. These strip graphs can be embedded on surfaces with the follow- 
ing topologies: (i) (FBC^FBOJ: strip; (ii) (PBC^FBO,): cylindrical; (iii) (FBC^PBQJ: 
cylindrical (denoted cyclic here); (iv) (FBC^TPBQJ: Mdbius; (v) (PBC^PBQJ: torus; 
and (vi) (PBC y ,TPBC x ): Klein bottle^] We have calculated P, W, and B for a variety of 
strip graphs having BC's of type (i) f30|- |3c 



(ii) |3[, (iii) H- f37j, (iv) |35|,|39|. Recently, 
strip graphs of the square lattice with BC's of torus (v) and Klein bottle (vi) type have also 
been studied Ifjlj. In addition to strips of regular lattices, we have calculated P, W, and B 



for cyclic graphs of polygons linked in various manners |37j] and homeomorphic expansions 
of strip graphs [[H].[||.[3(| . 

We proceed to describe our generating function method for calculating chromatic poly- 
nomials. The generating function is denoted. r(G s , q, x), and the chromatic polynomials for 
the strip of length L x = m are determined as the coefficients in a Taylor series expansion of 
this generating function in an auxiliary variable x about x = 0: 



T(G s ,q,x)= J2 P{(G s )m,q)x 



m—mo 



(2.1) 



m=mo 



where mo depends on the type of strip graph G s and is naturally chosen as the minimal value 
of m for which the graph is well defined. The generating functions T(G s ,q, x) are rational 
functions of the form p0| 

Af(G s ,q,x) 



with 



r(G s ,q,x) 



M{G S} q,x) 



V(G s ,q,x) 



(q)x J 



(2.2) 



(2.3) 



5=0 



These BC's can all be implemented in a manner that is uniform in the length L x ; the case (vii) 
(TPBCj,,TPBC x ) with the topology of the projective plane requires different identifications as L x varies 
and is not considered here. 
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and 



dp 
3=1 



.r 



(2.4) 



where the Ac St i and bc s ,i are polynomials in q (with no common factors). Writing the 
denominator of T(G S , q, x) in factorized form, we have 



V(G s ,q,x) 



deg x V 
3=1 



(2.5) 



One can then calculate Cj and aj in ( |1.4j ) in terms of these quantities: 

1 



djj r dj\f 

3=1 L s=0 



s-1 



n 



(Aj - A,) 



Aj 



(2.6) 



We have proved that for strip graphs with (FBC^FBCy) boundary conditions, all of the 
Aj's contribute to P, i.e. none of the coefficients in ( |2.6|) vanish [pO| , pT| . However, for some 
graphs with (FBC^PBCy) or (FBC^TPBC,,) boundary conditions, we have also proved that 
certain coefficients Cj in eq. (|2.6|) do vanish, so that the corresponding Xj;s do not contribute 
to P and iV a < deg x {V) [p^ , P9| . In general, for the A/s that do contribute, 



A, 



(2.7) 



We show here how the nonpolynomial algebraic roots in the various P's yield polynomials 
in q. As is evident from (|2.5|) , for a given strip, the A/s arise as roots of the equation T> = 0. 
In general, T> contains some number of factors of linear, quadratic, cubic, etc. order in x. 
Consider a generic factor in T> of r'th degree in x: (1 + fix + f'2X 2 + ... + f r x r ), where the fj's 
are polynomials in x. This yields r A^'s as roots of the equation £ r + /i^ r_1 + ... + f r — 0. 
The expressions in P involving these roots are symmetric polynomial functions of the roots, 
namely terms of the form 



{21 



and, for the Mobius strips of the L y = 2 triangular lattice and asymmetric homeomorphic 
expansions of the L y = 2 Mobius square strip, as well as strips with (FBC y ,FBC x ), terms of 
the form 



(A+r - (A-y 



m—1 



A, 



A_ 



^(A + ) m - x - fe (A_) fc . 



(2.9) 



k=0 
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We can then apply a standard theorem (e.g. ||44j| ) which states that a symmetric polynomial 
function of the roots of an algebraic equation is a polynomial in the coefficients, here fi, 
I = 1, ..,r; hence this function is a polynomial in q. For example, for s r , one has the well- 
known formulas (due to Newton) si = —fx, S2 = fx ~ 2/ 2 , S3 = —ff + 3/1/2 — 3/3, etc. 
It becomes progressively more and more time-consuming to calculate these s r 's for large m 
and hence large r; it is here that one uses the full power of the generating function method, 
which immediately yields the chromatic polynomial without the necessity of having to go 
through the intermediate stage of calculating the Ay's and then use Newton identities to get 
rid of algebraic roots and obtain the final polynomial. 

III. GENERAL RESULTS ON B 

One can ask a number of questions about the properties of B. We take the opportunity 
here to give a unified discussion of these which includes results from our various studies 
together with some new remarks. 

1. Does B have general symmetries? The answer is yes; a basic symmetry is that 

B(q) = B(q*) (3.1) 

i.e., B is invariant under complex conjugation in the q plane. This follows from the 
fact that the chromatic zeros have this property, which, in turn, follows from the fact 
that the coefficients in the chromatic polynomial are real (cf. eq. (|1.3|) ). 



What is the dimensionality of B in the q plane? As noted after (|1.4f) , one proves from 
the definition of B as the solution of the equality in magnitude of leading terms cij 
in eq. (|1.4j ), together with the fact that the a/s are algebraic functions of q, that (in 
cases where it is not the empty set) that B is an algebraic curve (including possible 
line segments on the real q axis), so 

dim(B) = 1 (3.2) 

What is the topology of the curve B and the associated "region diagram" in the q 
plane? Does B separate the q plane into two or more regions or not? Does B cross the 
real q axis, and, if so, at which points? From our studies [^0| , ^l| ,^5|- ^7l|39| , we arrive at 
the following inference, which we state as a conjecture: 
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Conjecture. For a family of graphs G s with well-defined lattice structure^], a sufficient 
condition for B to separate the q plane into different regions is that G s contains at 
least one global circuit, defined as a route following a lattice direction which has the 
topology of S 1 and a length £ ff . c . that goes to infinity as n — > oo. In the context 
of strip graphs, this is equivalent to having PBC X , i.e., periodic boundary conditions 
(or TPBCz, twisted periodic boundary conditions) in the direction in which the strip 
length goes to infinity as m — > oo. 

The presence of a global circuit in a family of graphs is not a necessary condition for B 
to enclose regions, as was shown in |32| . We concentrate here on lattice strips because 
of the connection to statistical mechanics; however, we also have studied families of 
graphs which do not have a lattice structure but, in the n — > oo limit, yield loci B 
that separate the q plane into regions. These include cyclic chains of polygons p5| , |37 



(for which the analogue of the global circuit is the route around the chain) and certain 



families of graphs with noncompact B p3|j2qj29|j3^j38|1 ; these are discussed below. A 
related conjecture is that, for a family of graphs G s with well-defined lattice structure, 
a necessary and sufficient condition for B to separate the q plane into regions and pass 
through q = is that G s contains at least one global circuit. Indeed, all of the families 
of graphs that we have studied that contain global circuits also yield loci B that pass 
through q = 2 as well as q = 0. In contrast, for strip graphs with (FBC^FBQz) studied 
in ||30| , |32"|1 , B consists of arcs that, in the simplest cases, do not enclose regions (c.f. 
Figs. 3-9 of [j30fl , Figs. 2(a), 3(a) of |32| , and in other cases do enclose regions (c.f. Figs. 



2(b), 3(b) of f32|), but do not pass through q = 0. We find that as the width of the strip 
increases, these arcs tend to elongate and move toward each other, thereby suggesting 
that if one considered the sequence of strip graphs of this type with width L y and 
length L x (having taken the limit L x — > oo to obtain a locus B for each member of 
this sequence), then in the limit L y — > oo, the arcs would close to form a B that passes 
through q = in such a way as always to separate the q plane into different regions. 
The interesting feature of the cyclic and Mobius strips and cyclic polygon chain graphs 
p5|- |37| , |39l is that when the graphs contain a global circuit, this property of B that it 



passes through q = in such a manner as to separate the q plane into regions already 
occurs for finite L y . This means that the W functions of these graphs with cyclic 
and twisted cyclic longitudinal boundary conditions already exhibit a feature which is 



By well-defined lattice structure we mean that the vertices and edges of the graph can be put into a 1-1 
correspondence with the vertices and edges of a section of a regular lattice. 
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expected to occur for the B for the W function of the full two-dimensional lattice. This 
expectation is supported by the calculation of W and B for the 2D triangular lattice 
with cylindrical boundary conditions by Baxter [ 24] . Recently, in |4l| with Biggs, we 



have found the same property for strips of the square lattice with (PBCj^PBCj;) (torus) 
and (PBC^TPBQr) (Klein bottle) boundary conditions. 

4. Does B have singular points in the technical terminology of algebraic geometry, such 
as (i) multiple points, where several branches of the curve intersect without crossing or 
cross each other, or (ii) endpoints? All of these possibilities are realized. There are cases 
where there are no multiple points, such as for the circuit and p-wheel graphs (defined 
in eq. (|3.3|) below) and some polygon chain graphs. There are also cases where there 
are (a) multiple points associated with intersecting but non-crossing branches, such as 
the L y = 2 and L y = 3 cyclic and Mobius strips of the square lattice and the L y = 2 
strips of the kagome (3 • 6 • 3 ■ 6) lattice P, |36| , |39H ; (b) multiple points associated with 



crossing curves, such as the L y = 2 cyclic and Mobius strips of the triangular lattices 
and certain families of cyclic polygon graphs; and (c) multiple points associated with 
both branch intersections and crossings, such for as homeomorphic expansions of cyclic 
square strips |35| . We have also shown that homogeneous strips of regular lattices with 



free transverse and longitudinal boundary conditions have loci B with endpoints |30 
(see further below). 

5. Does B consist of a single connected component, or several distinct components? In 
cases where B does not contain any multiple points, the number of regions, N reg . and 
the number of connected components on B satisfy the relation N reg , = N comp , + 1. The 



Harnack theorem ]45] gives the upper bound N comp . < h + 1, where h, the genus of the 
algebraic curve comprising B, is h = (d— l)(d — 2)/2. However, as we have discussed 
in 



34,37], this is a very weak bound. For example, for the cyclic polygon chain family 



denoted (ei,e 2 ,e 9 ) = (2,2,1) in [07H , we showed that N comp . = 2, while the Harnack 



theorem gives the upper bound N comp . < 37. 

6. Although B is the continuous accumulation set of a subset of the chromatic zeros, these 
zeros do not, in general, lie precisely on this asymptotic locus B for finite n. Are there 
families of graphs for which all or a subset of the zeros do lie exactly on B for finite 
n? We have answered this in the affirmative and have constructed a family of graphs, 
which we call p- wheels, and proved that except for a finite subset of chromatic zeros at 
certain nonnegative integer values (see below), all of the chromatic zeros lie exactly on 



B [ 25| . Recall that the join of two graphs G and H is defined as the graph consisting of 
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copies of G and H with additional edges added joining each vertex of G to each vertex 
of H ||. A p-wheel is the join 

(Wh)M = K P + C n . v (3.3) 

so that (Wh)n^ is the circuit graph, C n , (Wh)tP is the usual wheel graph, etc. The 
locus B is the circle 

|g-p-l| = l (3.4) 

and the real chromatic zeros occur at q — 0, 1, ...p+1 for n—p even and q = 0,1, ...,p+2 
for n — p odd, while the complex zeros all lie on the above unit circle. 

7. What is the density of the zeros on Bl In special cases, in particular, the p-wheel 
graphs, it is a constant. In general, it varies as one moves along B for a given family of 
graphs. Our calculations of chromatic zeros answer this question for specific families. 

8. For many years, no examples of chromatic zeros were found with negative real parts, 
leading to the conjecture that Re(q) > for any chromatic zero |[21| . Although Read 
and Royle later showed that this conjecture is false |23], very few cases of graphs with 
chromatic zeros having Re(q) < are known, and the investigation of such cases is 
thus valuable for the insight it yields into properties of chromatic zeros. Note that the 
condition that a graph has some chromatic zeros with Re(q) < is a necessary but 
not sufficient condition that it has an accumulation set B with support for Re(q) < 0. 
We have proved that certain homeomorphic expansions of square strip graphs with 



the boundary conditions (FBC y ,FBC x ) |3T[] and (FBC^PBCy |p5| , of length m units 



have chromatic zeros with Re(q) < for sufficiently large m and that they have loci B 



with support for Re(q) < 0. Since the homeomorphically expanded lattice strips in |3T 
do not have global circuits, this shows that these circuits are not a necessary condition 
for a recursive family to have chromatic zeros and B with support for Re(q) < 0, as 



is also shown for lattice strips with (PBC y , FBC X ) in [2J j, falsifying a conjecture in 
^C[. We have also shown the existence of chromatic zeros and B with Re(q) < for 
the L y = 3 cyclic and Mobius strips of the square lattice and the L y = 2 cyclic and 



Mobius strips of the kagome (3 • 6 ■ 3 ■ 6) lattice |36|j39|| . A different example is families 
of graphs with noncompact B [p|, p8| , p9| , p4| ; with these, we constructed cases where the 
chromatic zeros have arbitrarily large negative Re(q). 

9. Is B compact (which in our context is synonymous with the property of boundedness) 
in the q plane, or does it extend to complex infinity in some directions? This is related 
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to the question of boundedness of the magnitudes of chromatic zeros. In p3| Read 
and Royle gave an example of a graph with a noncompact B, namely the bipyramid 
graph. We have studied the question of the compactness of B in a series of papers 
[0.^,^,0 (see also |j38|). A necessary and sufficient condition such that B is compact 
in the q plane is obtained as follows: one simply reexpresses the degeneracy equation 
for leading a.,(g)'s in terms of the variable z = 1/q and determines if this has a solution 
for z — 0. By constructing and studying families of graphs with noncompact loci B, 
we have elucidated the geometrical conditions leading to this noncompactness. Our 
studies show that a necessary condition is that some vertex has a degree A which 
goes to infinity as n — ► oo. However, this is obviously not a sufficient condition for 
noncompactness of B, as is shown, for example, by the p-wheel graphs (Wh)^) n , in 
which the p vertices in the complete graph K p have degree A that goes to infinity as 
n — > oo but for which B is compact (eq. ( |3.4|) . From our studies we have found that 
in all cases of families with noncompact B, the graphs have the common feature that, 
in the limit as n — > oo, they contain an infinite number of different, non-overlapping 
non-self- intersecting circuits, each of which passes through two or more nonadjacent 
vertices. We are led to propose this as a conjecture for the condition on a graph family 



such that it has a noncompact B [|34]j2 



10. The above conjecture also leads to the following corollary: A sufficient (not necessary) 
condition for a family of graphs to have a compact, bounded locus B is that it is a 
regular lattice [p8[, ^9| , j34ll . Clearly, if B is noncompact, passing through the origin of 
the z plane, where z — 1/q, then the function q~ l W has no large-g expansion, i.e., no 
Taylor series expansion in the variable z around the point z = 0. Our conjecture is in 



accord with the derivation of large-g expansions for regular lattices ]T7| 



11. Another feature that we find is that for families of graphs that (a) contain global 
circuits, (b) cannot be written as the join G = K p + H, where K p is the complete 
graph on p vertices, and (c) have compact B, this locus passes through q = and 
crosses the positive real axis, thereby always defining a q c . Note that B({K P + H}, q) = 
B({H}, q — p), i.e. the B for {K p + H} is the same as that for {H} shifted a distance 
p to the right in the q plane. Hence, examples of families of graphs with loci B that 
do not pass through q = include these joins, where if is a family whose B does pass 
through zero. Other examples of graph families with B not including q = are those 
with noncompact B |p8| , |29| , |34l . 

12. What is the effect of the boundary conditions of the strip graph G s on £>? For homo- 
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geneous strip graphs with (FBC^FBCy boundary conditions, in the cases studied in 
[30H where B is nontrivial, it consists of arcs and hence does contain endpoint singular- 



ities. In these cases, the X/s are algebraic expressions and the arcs comprising B run 
between branch point zeros of algebraic roots of polynomials in q. An arc can cross 
the real axis q axis if and only if it is self-conjugate. B contains a line segment on the 
positive real axis, say in the interval q± < q < q2, if there are two X/s, e.g., of the form 
\jj> = R(q) ± a/ S(q) with R and S being polynomials in q such that S < in the 
interval q\ < q < q 2 . 

13. Again concerning the question of the role of boundary conditions: for the strip graph 
{Gs)m with a given type of transverse boundary conditions BC^, the chromatic poly- 
nomial for PBCz has a larger N a than the chromatic polynomial for FBC X , and the 
corresponding loci B are different. 

14. In the examples we have studied, we find that for a given type of strip graph G s with 
FBC y , the chromatic polynomials for PBC X and TPBC X boundary conditions (i.e., 
cyclic and Mobius strips) have the same aj, although in general different Cj. It follows 
that the loci B are the same for these two different longitudinal boundary condition 



choices ]35]|36]|39 



15. However, in the case of PBC y , the reversal of orientation involved in going from PBC a . 
to TPBCz longitudinal boundary conditions (i.e. from torus to Klein bottle topol- 
ogy) can lead to the removal of some of the a^s that were present; i.e., P for the 
(PBCjpTPBC-j;) strip may involve only a subset of the a/s that are present for the 
(PBC^PBC-j;) strip. For example, for the L y = 3 strips of the square lattice with 
(PBCy, PBC X ) boundary conditions, there are N a = 8 Oj's, but for the strip with 
(PBCy, TPBC X ) boundary conditions only a subset of N a = 5 of these terms occurs in 
P [f|T|. None of the three a/s that are absent from P in the TPBC^ case is leading, 
so that B is the same for both of these families. 

16. How do W functions in region Ri compare for different boundary conditions? We have 
shown that, for a given type of strip graph G s , in the region R\ defined for the PBC^ 
boundary conditions, the W function is the same for FBC^, PBC X , and TPBC^. 

17. For strips of regular lattices, as one increases L y , how does W approach the limit for the 
2D lattice? We have found that with q > q c for a given lattice type A, the approach of 
W to its 2D thermodynamic limit as L y increases is quite rapid; for moderate values of 
q and L y ~ 4, W(A, (L x = oo) x L y , q) is within about 5 % and C(10 -3 ) of the 2D value 
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W(A, [L x = oo) x (L y = oo), q) for FBC y and PBC y , respectively for a strip graph of 
the lattice A. We have proved that the approach of W to the 2D thermodynamic limit 
is monotonic for FBC^ and non-monotonic, although more rapid, for PBC y ||33|| . (By 
the result noted in the previous item, the W function for these values of q, in region 
Ri, is independent of the longitudinal boundary conditions.) 

18. Concerning the analytic properties of W, one property is that W may have isolated 
branch point singularities (zeros) for strips with (FBC^FBCj;) and (PBC y ,FBC x ). For 
example, for the (FBC y ,FBC x ) strips of the square lattice with L y = 2, and L y = 3, 
the respective W functions have square and cube root branch point singularities, while 
for the (PBC^FBCa;) strips of the square lattice with L y = 3, i.e. transverse cross 
sections forming triangles, the W functions have cube root branch point singularities 
where they vanish. Note that for these cases, Ri is the full q plane, and the continuous 
locus B — 0. In contrast, for the corresponding lattice strips PBC^, although the W 
functions have the same expressions in the respective R\ regions, given below as eqs. 
1|), ( f4.10|) , and (|4.22|) , the zeros of the expressions do not occur in the R\ regions 



where these expressions apply, and so the W functions in these R\ regions do not have 
any isolated branch point singularities. We have also found this to be true of other 
strips. 

IV. CHROMATIC POLYNOMIALS AND LOCI B FOR SOME SPECIFIC 

FAMILIES 

Here we present exact calculations of chromatic polynomials, W functions, and loci B for 
various families of graphs. 

A. Strips with Free Transverse and Longitudinal Boundary Conditions 

We first show B for the strip graphs of the square lattice with (FBC^FBC^) and L y = 2, 3. 
The relevant generating functions are given in ||3(| . 
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FIG. 1. B for W{{G sq ( Ly )},q) with L y = (a) 3, (b) 4, where {G sq ( Ly )} denotes the (L x = oo) X L y 
strip of the square lattice. For comparison, the zeros of the chromatic polynomial P({G sq (i, )) m , q) for (a) 
L y = 3, m = 16 (hence n = 54 vertices) and (b) L y = 4, m = 16 (hence n = 72) are shown. 

B. Strips with Free Transverse and Periodic or Mobius Longitudinal Boundary 

Conditions 

For definiteriess, we consider cyclic and Mobius strips of the square lattice. For L y = 1, 
these reduce to the circuit graph, for which B is the circle \q — 1| = 1. For the L y = 2 strips 
of this type, from the P in || we obtain B given in Fig. |[ 
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FIG. 2. B for the infinite-length limit of strips of the square lattice with (FBCj^PBCz) (cyclic) boundary 
conditions and L y = (a) 2, (b) 3. For comparison, the zeros of the chromatic polynomial P((G S9( x y )) m , q) 
for (a) to = 19 (hence n = 38) and (b) m = 20 (hence n = 60) are shown for comparison. The B's for the 
corresponding strips with (FBC y ,TPBC x ) (Mobius) boundary conditions are identical. 



For the cyclic L y = 3 strip of the square lattice, we find j3(^ 
P(sq(L y = 3, cyc.) m , q) = (q 3 - 5g 2 + Qq - l)(-l) m 



+tf - 3 q + 1) ( g _ 1)- + (g _ 2 ) m + (g - 4) m + (g - l)[-(g - 2) 2 ] 



\2im 



+ 



(A sg , 6 ) m + (A S9i7 r +(9-1) (A S9l8 ) m + (A^ 9 ) m + (A s ,,io) 



where 



As<?,(6,7) 



(g - 2)(g 2 - 3g + 5) ± <{ (g 2 - 5g + 7)(g 4 - 5g 3 + llq 2 - 12q + 



1/2' 



(4.1) 



(4.2) 
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and X S q,j, j = 8, 9, 10, are the roots of the cubic equation 



with 



b sq ,2i = 2q 2 - 9q + 12 



b sq ,22 = q 4 ~ 10g 3 + 36q 2 - 56q + 31 



(4.3) 

(4.4) 
(4.5) 



b SQ ,23 = ~{q ~ l)(g 4 - 9g 3 + 29g 2 - 40g + 22) . (4.6) 
and for the corresponding Mobius strip |39| 

P(sq(L y = 3, Mb.) m , q) = (q 2 - 3q + l)(-l) m - (q - l) m + (q - 2) m - (q - 4) m 
-(9 " l)h(? " 2) 2 ] m + [(X sq ,6) m + (A sgJ ) m ] + (q - 1) [(A sg , 8 ) m + (X sq , 9 ) m + (X sqAO y j 

(4.7) 

For the (L^ — > 00 limits of the) L y — 1 and L y = 2 cyclic strips, q c = 2, while for 
the L y = 3 cyclic strip, g c = 2.33654. These values are obtained from exact solutions for 
the respective W and B for these strips P,PBL3D1. For these cyclic graphs this point q c is a 
non- decreasing function of L y . In the limit L y —* 00, continuity arguments imply that q c 
approaches the value for the 2D square lattice, which is q c (sq) = 3 [ |i~Tf . 

The W functions for the regions are given in [B,5B|,RD|. In particular, note that in the 
respective region R\ for each family, 



W{sq(L y = 1), eye, q) = q - 1 
and, for both cyclic and Mobius strips, 

W(sq(L y = 2), eye, q) = (q 2 - 3q + 3) 1 / 2 

W(sq(L y = 3, eye), q) = 2" 1 / 3 [(g - 2){q 2 - 3q + 5) + 

(q 2 -5q + 7)(q 4 - 5q 3 + llg 2 - I2q + 



(4.8) 



(4.9) 



1/2 



1/3 



(4.10) 



These W functions are the same as for the corresponding strips with (FBC y ,PBC x ), which 
is a general result, as noted above. 
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C. Homeomorphic Expansion of Cyclic and Mobius Strips 



For comparison with the previous figure, we show homeomorphic expansions obtained by 
starting with cyclic and Mobius L y = 2 square strips consisting of m squares, and inserting 
k — 2 additional vertices on each horizontal edge. We denote these graphs as as {Ch) k , mfiyc . 
and (Ch) k,m,Mb.- These graphs can be regarded as cyclic and Mobius strips of m p-sided 
polygons, where p = 2k, such that each p-gon intersects the previous one on one of its edges, 
and intersects the next one on its opposite edge. It is convenient to define 



fc-2 



where P(C m , q) is the chromatic polynomial for the circuit (cyclic) graph C m with m vertices, 
P(C m , q) = {q~ l) m + {q- l)(-l) m We calculate |5| 

P{{Ch) k , m ,cyc, q) = q 2 -3q + l + (D 2k ) m + 

- / \ m / \ in- 

(q-1) ((-l) k+1 D k+1 + D k ) +((-l)^D k+1 -D k ) (4.12) 



P((Ch) 



k,m,cyc.,Mb 



,q) = -l + (D 2k ) m + 



-lY(q-l) [(-l) k+1 D k+1 + D 



\k+l 



D 



fe+i 



Di 



(4.13) 



In Fig. ^ we show B and chromatic zeros for the cases k = 3,4. In region Ri, which includes 
the real axis for q > 2, 



W = (D 2k )w=V 



(4.14) 
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D. Cyclic Chains of Polygons 



Consider a cyclic chain composed of m subunits, each subunit consisting of a p-sided 
polygon with one of its vertices attached to a line segment of length e g edges (bonds). Thus, 
the members of each successive pair of polygons are separated from each other by a distance 
(gap) of e g bonds along these line segments, with e g = representing the case of contiguous 
polygons. Since each polygon is connected to the rest of the chain at two vertices (taken to 
be at the same relative positions on the polygons in all cases), this family of graphs depends 
on two additional parameters, namely the number of edges of the polygons between these 
two connection vertices, moving in opposite directions along the polygon, e% and e<i- We 
denote this family of cyclic polygon chain graphs as G eije2 ,e B ,m;a and define, in this context, 
p = e\ + 62- A symmetry property is 



We calculate |37] 



P(Ge lt e 2 ,e g ,m: <Z) = K)™ + (<7 ~ !){<h) m (4.16) 

where 

a l = {q-iy^ l D p (4.17) 



a 2 



-i 



q-2 + (1 -g) ei + (1 



■1) 



p+e g 



ei 

£ 

s=2 



q) 

E 



-qY 



(4.18) 
(4.19) 



The locus B obtained by taking the limit m — > oo is shown in Fig. ^ for the cases 
(ei,e 2 ,e 9 ) = (a) (2,2,0), (b) (2,2,1). One sees the interesting phenomenon that as the 
number of edges in the gap e g increases from 0, the multiple point that was present on B for 
e g = disappears and B decomposes into two separate components. As e g increases further, 
the inner boundary shrinks monotonically around the point at its center. This family also 
illustrates the fact that one can take the limit n — > oo in a different way, letting ei, e2, or 
e g go to infinity with m held fixed. For the nontrivial case where min{e\,e<2) > 1, we find 
that, if p is even, then q c = 2 while if p is odd, then q c < 2 and, for fixed (ei, e 2 ), q c increases 
monotonically as e g increases, approaching 2 from below as e g — > oo. In region R 1: 



W = [{q-iy^Dp}^- 1 



(4.20) 
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which is again the same W function as for the corresponding polygon chain graph with free 
(open) longitudinal boundary conditions. 
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FIG. 4. Boundary B in the q plane for W function for linim^oo G ei ,e 2 ,e g ,m with (ei, e2, e g ) = (a) (2,2,0), 
(b) (2,2,1). Chromatic zeros for m = 14 (i.e., n = 42 and n = 56 for (a) and (b)) are shown for comparison. 

E. Families with Periodic Transverse and Free Longitudinal Boundary Conditions 

Here we consider strip graphs with (PBC^FBCy. For the strip graph of the square 



lattice with L y = 3, i.e., cross sections forming triangles, and FBC X , we find ||30| , |33 

P(sq(L y = 3) m , PBC y , FBC X , q) = q(q - l)(q - 2)(q 3 - Qq 2 + Uq - 13)" 1 - 1 (4.21) 
whence 

W(sq(L y = 3), PBC y , FBC X , q) = (q 3 - Qq 2 + Uq - 13) 1/3 (4.22) 
with B — 0. Results for this lattice strip with larger L y = 4, i.e. cross sections forming 



squares, and for the triangular lattice, are given in [33]. 



F. Families with Periodic Transverse and (Twisted) Periodic Longitudinal Boundary 

Conditions 

With N. Biggs, we have recently studied the case of the strip of the square lattice with 



(PBC^PBCa) (torus) and (PBC^TPBC^) (Klein bottle) boundary conditions; see Ref. gl 
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for the chromatic polynomials and W functions. We find q c = 3, which, interestingly, is the 
same as for the full 2D square lattice. 



G. Families with Noncompact B 



As an example of a family of graphs with a noncompact B is obtained as follows [2S|. 
We start with the join (K p ) + G r and remove some bonds from K p . The simplest case is to 
let G = K r , K p = K 2 , and remove the edge connecting the two vertices of the K 2 . We can 
homeomorphically expand this by adding degree-2 vertices to the r edges connecting each 
of the two vertices of the original K 2 to the r vertices of the K r . We denote this family as 

§§. we find mm 



H, 



k.r 



P(H k , r ,q)=q(q-V 



D 



2fe-2 



\r-2 



r-2 



q(q-l) {q-iy-\D k ^y+{D 



(4.23) 



Because of the noncompactness of B in the q plane, it is more convenient to plot it in the 
- 1/q plane. In Fig. || we show two typical plots. For further details, see Refs. [p8jp9j^4 



General bounds o-fi chromatic zeros of graph have been discussed in a number of papers, e.g. 
and recently @. 
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FIG. 5. Boundary B in the z = l/q plane for lim r _>oo Hk, r with k = (a) 4 (b) 5. Chromatic zeros for 
Hk tr with (k,r) = (a) (4,30) (b) (5,18) are shown for comparison. 

V. ZEROS AND THEIR ACCUMULATION SETS FOR POTTS PARTITION 
FUNCTIONS AT FINITE TEMPERATURE OR TUTTE DICHROMATIC 

POLYNOMIALS 

We have recently generalized our studies to the case of the g-state Potts model at arbitrary 
temperatures. At temperature T on a graph G or lattice A this model is defined by the 
partition function 



where Oi = l,...,q are the spin variables on each vertex i e G; (3 = (ksT) 1 ; and (ij) 
denotes pairs of adjacent vertices. (Our terminology for this partition function is standard 




(5.1) 



with the Hamiltonian 




(5.2) 
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in physics; obviously Z should not be confused with the partition function p(n) in number 
theory and combinatorics.) We use the notation K = j3J, 



a = u~ l = e K (5.3) 

and 

v = a-l (5.4) 

and denote the (reduced) free energy per vertex (site) as / = — j3F = lim^oo rT 1 In Z. The 
partition function can be written as 

Z(G,q,a) = J2H( 1 + v5 ^ ( 5 - 5 ) 

R} (v) 



which shows that it is a polynomial in q and v or a. From ( |5.5| ), it follows that [ p^ , p^ ,p^,[7| 

Z(G,q,a)= q HG,) v e ^ (5.6) 
G'CG 

where G' is a subgraph of G and e{G') and k{G') denote the edges (bonds) and connected 
components, including single vertices, of G', respectively. 

The ferromagnetic and antiferromagnetic signs of the spin-spin exchange coupling are 
J > and J < 0, respectively, and hence as the temperature varies from to oo, the 
variable a varies from to 1 for the Potts antiferromagnet and from oo to 1 for the Potts 
ferromagnet. At T = oo, i.e., /3 — or a — 1, the Potts ferromagnet and antiferromagnet 
are identical, since J does not enter in Z, which reduces simply to 

Z(G,q,a = l) = q n ^ (5.7) 

As indicated in eq. (|1.1|) , the T = (a = 0) case yields the chromatic polynomial. For other 
special cases one has the elementary results Z(G, q = 0, a) = and 

Z(G,q = I, a) = a e{G) (5.8) 

The zeros of Z are now functions of both complex q and a. In the limit n —* oo, these zeros 
have an accumulation set which is a singular submanifold B in the C 2 space defined by q 
and a (or equivalently, it). Just as, in the one- variable situation, the W({G}, q) function was 
singular across the locus B in the q plane as it switched its analytic form due to a change in 
the dominant term a,j(q) in (|1.4j ), so in the two- variable situation, the reduced free energy is 
singular across the locus B in the (q, a) space as it switches it analytic form due to a change 
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in the dominant term Xj in (|5.13|) . We denote the slices of this submanifold B in the q plane, 



and in the a or u variable as B q , B a , and B u . 



We recall that the partition function Z(G, q, a) is related to the Tutte (dichromatic) 
polynomial T(G,x,y) |[4|,[yJ and the rank function R(G,£,,rj) fl3|| . Defining 



x=l + - (5.9) 

v 

and 

y = a = v + l (5.10) 

so that q = (x — l)(y — 1), one has 

T(G, x, y) = (x- iy k{G \y - l)' n ^Z(G, q, v) (5.11) 

This is equivalent to the usual expression of the Tutte polynomial in terms of spanning trees 
||14j| . The connection with the rank function follows from the relation T(G,x + l,y + 1) = 
x n ~ 1 R(G, x" 1 , y), viz., 

Z(G, q, v) = q n ^R(G, £ = -, V = v) (5.12) 

Q 

The partition function or Tutte polynomial can be calculated either by the iterative use of 
the deletion-contraction theorem or by a generalization of the usual transfer matrix method 
from statistical mechanics. For recursive families, from the transfer matrix method, it follows 
that Z has the structure 

Z(G,q,a) = J2cAXjr (5-13) 

where the Cj and the Xj depend on G but are independent of m. This is a generalization 
of eq. ( |1.4| ) (with the equivalent notation N\ = N a and Xj = aj). Given the formula ( |5.6| ), 
the zero-temperature limit of the Potts antiferromagnet is studied by taking a — >• 0. For 
the Potts ferromagnet, since a — > oo as T — > and Z diverges like a e ( G ) in this limit, it is 
convenient to use the low-temperature variable u = 1/a = e~ K and the reduced partition 
function Z r defined by 

Z r = a~ e{G) Z = u e{G) Z (5.14) 

which has the finite limit Z r — > 1 as T — > 0. For a general strip graph (G s ) m of type G s and 
length L x = m, we can write 
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Z r ((G s ) m ,q,a) = u<^ J>(A;r = E c ^'(V) m (5-15) 



3=1 3=1 

with 



\ jr = M e « Gs ^)/ m A, (5.16) 



An elementary calculation for the circuit graph C n (e.g., |46[]) yields Z[G = C n ,q,v) = 
(q + v) n + (q — l)v n . The locus B is the solution of the degeneracy equation \a + q — 1 1 = \a — 1| 
£>. The slice in the q plane, B q , thus consists of the circle centered at q = 1 — a with radius 
|1 - a\: 

q= (l-a)(l + e ie ) , 9e[0,2n) (5.17) 

For the Potts antiferromagnet, as a increases from to 1, this circle contracts monotonically 
toward the origin, and at a = 1 it degenerates into a point at q — 0. For the Potts 
ferromagnet, as a increases from 1 to oo, the circle expands into the Re(q) < half-plane. 
One can also consider other values of a that do not correspond to physical temperature in 
the Potts ferromagnet or antiferromagnet. In all cases, from eq. ( |5.17| ) it is evident that B q 
passes through the origin, q = 0. For real a ^ 1, B intersects the real q axis at q = and at 

q c ({C}) = 2(1 - a) (5.18) 

If, as is customary in physics, one restricts to q G Z + , then the Potts antiferromagnet on {C} 
has a zero-temperature phase transition if and only if q — 2 (Ising case). In the ferromagnetic 
case, B q does not cross the positive real q axis. For q ^ 0, 2, the slice of this submanifold in 
the u plane forms the circle 

B u : u= (q-2)- 1 (-l +0 , cjG[0,2tt), q ^ 0, 2 (5.19) 

while for q = 2, B u is the imaginary u axis. For q = 0, Z = and no B u is defined. We 
have obtained similar results for other graphs, including cyclic ladder graphs L m ||42j| . The 
corresponding Tutte polynomial obeys a recursion relation given in |J. Some zeros in the q 
plane are shown in Figs. |] and [7]: 
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FIG. 6. Zeros of Z{L mi q, a) in q for m = 18 and a = 0.25. The axis labels qr = Re(q) and qi = Im(q). 
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FIG. 7. Zeros of Z(L m , q, a) in q for m = 18 and a = 0.5. 

The locus B always passes through the origin, q = 0. For < a < 1, B crosses the 
positive real q axis at q c , where 

q c ({L}) = (l-a)(a + 2) for 0<a<l (5.20) 

and separates the q plane into several regions. As a increases from to 1, the locus B 
contracts monotonically toward the origin, q = and in the limit as a — > 1, it degenerates 
to a point at q — 0. This also describes the general behavior of the partition function 
(dichromatic) zeros themselves. That is, for finite graphs, there are no isolated partition 
function zeros whose moduli remains large as a — > 1. This is clear from continuity arguments 
in this limit, given eq. ( |5.7| ). 

VI. RIGOROUS BOUNDS ON W 

Using a coloring matrix method, Biggs [19|] proved upper and lower bounds for W for 
the square lattice. We applied this method to prove such bounds for the triangular and 
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honeycomb lattices [2(| . These bounds are very restrictive even for moderate q. Although a 



bound on a given function need not, a priori, coincide with a series expansion of that function, 
the lower bound for the honeycomb lattice coincides with the first eleven terms of the large-g 
expansion for W(hc, q). We have proved a general lower bound, which is applicable for any 
Archimedean lattice ||271| . Here an Archimedean lattice is defined as a uniform tiling of the 
plane with one or more regular polygons such that all vertices are equivalent to each other. 
It can be specified by the ordered sequence of polygons pi traversed by a circuit around any 
vertex: A = YliPi 1 ■ Let J2 a i = a M an d v% = a,i,s/Pi- Then our general lower bound is ]27| 



For the function 



this bound reads 



W > 



W(A,y) 



W(A,q) 
q{l-q- 1 )^/ 2 



w > n[ 1 + {- l ) n y n ~ 



(6.1) 



(6.2) 



(6.3) 
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